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SUMMARY 


Magnetohydrodynamics theory has been used to study the dynamic response of 
the neutral atmosphere to a geomagnetic storm. A full set of magnetohydrodynamic 
equations appropriate for the present problem is derived and their various orders of 
approximation are discussed in some detail . In order to demonstrate the usefulness 
of this theoretical model/ the May 1967 geomagnetic storm data recorded at College, 
Alaska; Dallas, Texas; and Honolulu, Hawaii have been used in the resulting set of 
non-linear, time dependent, partial differential magnetohydrodynamic equations to 
calculate variations of the thermosphere due to the storm. The numerical results are 
presented for wind speeds, electric field strength, and amount of joule heating at a 
constant altitude for the data recorded at the above-mentioned stations. They show 
that the strongest thermospheric responses are at the polar region (i »e». College, 
Alaska), becoming weaker in the equatorial region (i.e., Honolulu, Hawaii). This 
may lead to the speculation that a thermospheric wave is generated in the polar 
region due to the geomagnetic storm which propagates towards the equator. 
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CHAPTER I. 

INTRODUCTION 

A geomagnetic storm is a natural disturbance of the geomagnetic field on and 
above the earth's surface. Its origin and associated phenomena are very complicated. 
Theoretically, a variation of the magnetic field may be caused by the following rea- 
sons: 1) variation of the geomagentic poles, 2) presence of local nonhomogeneous 
electric field ■ ■■■ — - - V x E , and 3) variation of electric currents. 

O i 

The magnetic field generated by the earth's magnetic poles is relatively steady. 
Its variation is negligible in the period of a geomagnetic storm. On the surface of 
the earth the effect of the ionization of the atmosphere is small and there is no signifi- 
cant electrical field that can be detected. Therefore, it is believed that the geomag- 
netic disturbance recorded on the earth's surface is mainly due to the existence and 
variation of electric currents flowing above the earth (Chapman and Bartel [1 ] )• 
Indeed, if the electric currents data were given for the whole environment of the 
earth, then the magnetic field could be calculated everywhere, at least theoretically. 
This would explain the disturbance of the magnetic field recorded on the earth’s surface 
during the geomagnetic storm. Actually, there are always electric currents flowing 
above the earth, and their intensity and directions are changing from time to time. As 
a consequence, the geomagnetic field values recorded on the earth will never be 
constant. If the transient geomagnetic field variations are smooth and regular, we say 
it is quiet; otherwise it is described as disturbed. 

On a quiet day, the magnetic variation proceeds mainly according to local 
solar time with a small part related to the moon. The two parts are called solar daily 
and lunar daily magnetic variations, and their correspondence to magnetic fields is 
denoted by Sq and L. During a magnetic disturbance, additional electric currents 
flow in the ionosphere. They are superimposed on the Sq and L currents. From the 
analysis of geomagnetic disturbance data, we see that at least five components of 
electric currents are involved in a magnetic storm. They are described as follows: 
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1. DCF = disturbance due to current attributable to solar corpascular flux 
at magnetopause, 

2. DR = disturbance due to magnetic ring current, 

3. DP = polar current which is strongest in auroral electrojets, 

4. DT = disturbance due to magnetospheric tail currents, and 

5. DG = disturbance due to induced ground currents. 

Since the electric currents are flowing in the slightly ionized gas medium of 
the upper atmosphere, the other associated physical and dynamic effects such as wind 
generation, joule heating, and density and temperature variations can be predicted 
by magneto-gas-dynamic theory (Piddington [2 ]). A theoretical feature of the geo- 
magnetic storm and its associated phenomena are illustrated in Diagram I. (See also 
[2], P. 13.) 

Therefore, a complete analysis of the problem of the dynamic structure of the 
upper atmosphere due to a magnetic storm has to include consideration of the distur- 
bance currents and their associated electric and magnetic fields and the dynamic quali- 
ties simultaneously since they are interrelated to each other (see Diagram I). However, 
the solution of this problem is difficult due to the complexity of the mathematics. Many 
authors have studied one or a few particular effects separately and made assumptions 
about the other physical qualities rather arbitrarily; Cole 5] studied the joule heating 
of moving ionized gas on the assumption of steady and uniform electric and magnetic 

fields. He concluded that the joule heating effect is significant dueing a geomagnetic 

-5-3-1 

storm. The joule heating may be on the otder of 10 erg cm sec in the region of 

5 

100 to 200 km, and the wind as high as 10 cm/sec, if the assumed magnetic and elec- 
tric field configurations are correct. Later, he [4 ] extended the results by including 
viscous effects. Maeda and Kato [5 ]have given an excellent review on the problems 
of electrodynamics of the ionosphere in which the problems of conductivity, wind and 
the dynamo theory, drift and its effect on the ionospheric formation, and the interaction 
between wind and electromagnetic field are discussed in detail . Thomas [6] and 
Thomas and Ching [7 ], applying a one-dimensional vertical model, reproduced the 
height profile and the mean time lag of the density disturbance by assuring that the 
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heat input due to a magnetic storm is given. Volland and Mayer [8 ] re-analyzed the 
same problem using a three-dimensional thermospheric model . 

In this study, a theoretical model is established for calculating the joule heat- 
ing and winds from the geomagnetic variations recorded at storm time. The mathe- 
matical formulation is based on the magnetohydrodynamic theory (Chang, Wu, and 
Smith [9]). Faraday's law is employed for determining the electric field from the 
magnetic field disturbance data dueing the geomagnetic storm. 
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CHAPTER II. 

FORMULATION OF THE PROBLEM 
II— 1 . Basic Equations 

The thermosphere may be considered as a continuum medium with a finite 
electric conductivity (Cowling [10]). The presence of electric and magnetic fields 
in the conducting medium will give rise to two principal effects: First, body force 
(Lorentz force) and, second, energy generation (Joule heating). These must be taken 
into consideration in the momentum and energy equations. We will derive the equations 
of mass, momentum, and energy conservation as follows: 

Mass Conservation 

The equation of mass conservation is the same as in ordinary fluid dynamics, 

namely 

- + V • ( p v ) = 0 (1) 


where p denotes the mass density and v the velocity vector. 

Momentum Conservation 

The equation of motion of a continuum medium in general can be written 
as (Cauchy equation) 


P 


D v. 

dT 


ap. . 

LL 

B x . 
J 


+ F. 

i 


( 2 ) 


where v. are the components of the velocity vector v, P.. are components 
of the stress sensor, and F. denotes the components of the body force F. For a 
Newtonian fluid the stress tensor can be expressed as follows 


P.. 

•J 



V i 

V )6.. - T1 (- 

IJ 0 X. 

J 


Bv. 


+ 


Bx. 


J-) 


(3) 


where p is the static pressure and T] is the viscosity. By substituting Eq. ( 2 ) 
into Eq. ( 3 ) we obtain 

P =-Vp+x + F (4) 
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where 


X = ? ( T] 7 • v ) + T| [7 2 v + 7 (7 * v ) ] 


+ 2 C((v n) * v) v + (v n ) x (v x v) ] 


and the body force F is 


F - p g E + j x B + p g 

t — * 

where p^ is the charge density, E is the electric field strength, B is the 
magnetic induction, and j is the sum of the conduction current and the current 
flow due to convective transport of charges. 

Energy Conservation 

The rate of increase of total energy In the fluid of a moving volume, a , 
is given by 


D$ 
D t 


d a 


-L a +/. -2^ e) d a . 

2 D t { D t 


where e is the internal energy per unit mass. 

In order for energy to be conserved, this must equal the energy inputs per 
unit of time from other sources. These are 

1 . Joule heating = J (E • j ) d O’ 

2. Heat conduction per unit time = - J V* ( X. 7T )d a , where \ 

is the thermal conductivity. 

3. Work done by surface force = “ / 2L* v. P.. d s. applying Gauss’ 

J : r * 'J j 


theorem 


'/ J 


- f S v. P.. dS. ~ 

J i,j ' 'J J 


-/ ? 

i I 


77 Bx. "Vij 

' J J 


( v. P.. ) d a 


The energy equation is obtained as 


/ "5T d ° = /l-ido./v-f 
V \/ V 


X 7 T ) d a 


J v ( v. P. . ) d a 

V i/j 9x j 1 *J 
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( 5 ) 



The last summation on the right hand side of £q. (5) can be rewritten as 

- Xtt <p*i> + * 

* i 

where, from Eq. (3) 



Maxwell Equations and their Approximation 

The Maxwell equations are needed to determine the electromagnetic 
quantities. Thus, 


B = 0 


v • E 


v x E = 


SB 

St 


(6) 

(7) 

(8) 


-* -* ftp 

VxB=n(J+e-|£- 
o or 


(9) 

(10) 


and Ohm's Law; J-a(E+vx B ) + a 2 j^Bx (E + v x B)J/| B | 

A number of approximations can be made which will -be valid as long as the present 
analysis is concerned: 
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3 E 


1) The displacement current e - --— may be neglected in Eq«_(9). In order 

Of -v r 


3 P 

3 t 


to the conduc 


to show that is a valid approximation/ we shall compare e 
tion current J = cr E. If E is assumed to vary periodically with time with frequency uj 
then the ratio of the amplitude of e to the amplitude of J is approximately 

•(-m 


e ( 3 E / 3 t ) max _ e u) 
( a E ) max a 

.-3 


W 


/max 


max 


euu 


The value for o in our case is about 10 mhos/m. The value of e in a vacuum 

-12 -2 -1 

is approximately 9 x 10 farads/m. uu is about 10 sec . This gives 


e ( $ E / at) 


max 


= 10 


-10 


( o E )max 

which shows that the displacement current can be indeed neglected. 

2) The electrostatic body force 

If E varies linearly over some small region, then from Eq. (7) we see 


p = _iL 
e L 


Therefore 


r ~ 

P e E = 


eE 


In comparison to the Lorentzian force J x B = cr v B , we have 

P E _ 2 2 „2 

e ~ e E ~ evB _ ev 


JxB a Lv B 2 


aL v B 


2 a L 


The values in our problem are estimated as follows 

e = 10 11 farads/m 

3 / 

v - 10 m/sec 
L = 10^ m 


a - 10 mhos/i 


m 


v . 10 -11 


L a 

Therefore, the electrostatic body force can be neglected in comparison to the Lorentz 
force. 
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Now, let us recollect the above derived equations as follows: 
Maxwell 's equations: 


V x E = - 


a B 
a t 


( 11 ) 


V X B = V j 
o 


7 • B = 0 


( 12 ) 

(13) 


e V * E _ P 


(14) 


Conservation of Mass: 


Equation of Motion: 


- A P - + 7 

a t 


D v - 


(P v ) = 0 


D t 


- Vp + X + j x B 


(15) 


(16) 


Conservation of Energy: 

„ D $ 

P Dt 


E • j + V(\VT) -?*(pv)+<P 


(17) 


Ohm's Law: 


j - ", E' + °2 


B x E’ 

TbT 


(18) 


* — * “♦ 

where E ' = E+ v x B. This set of equations will form the basis of our theoret- 
ical analysis. Various approximations will be introduced as long as the physical 
situations permitt. 
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II— 2 . Some Fundamental Solutions 


Since the magnetohydrodynamic equations [Eq. (11) - Eq. (18) ] combine 
the full complexity of Maxwell's equations and the fluid dynamic equations, it is 
obvious that they will be extremely difficult to solve in the general form. Exact 
solutions exist only for a few special cases. We will discuss two of them; namely, 
the Piddington [ 11 ] solution for a non-viscous gas and the Hartman solution for a 
viscous fluid [12], They do not necessarily correspond to the conditions during mag- 
netic storm; however, the solutions will give us some insight into the physical nature 
of the problem. 

(A) Consider a uniform ionized gas in a steady and homogeneous electric 
field. They are assumed in the following forms 

E = E (E , E , 0 ) and B = (0, 0, B ) 
x y z 

We further assume that the pressure gradient is 


Vn = 


P " ' 


_iE_ 

3x » 


where i denotes the unit vector in the 


x-di recti on. 

By applying Eq. (18) and Eq. (16) and neglecting the viscous term, we 

obtain 

— ♦ — ♦ 

+ a ■ L 


j =a i E ‘ + °2 B 


(19) 


and 


o = 1 X 1 - 7p 


( 20 ) 


Dv 


At hydrostatic equilibrium = 0 > ^ en 


1 _ - V p X B _ 


_ , B x E’ 

- a E + o 

1 2 B 
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It follows 



2 

where cr^ = cr ^ + cr^ / • 

— * —♦ — ♦ — ♦ 

Therefore, the current j and velocity v are solved in terms of E, B and Vp . The 

Joule heating is given by 

Q = j* • E' = } • (E+ v x B) (23) 


This means that if we simulate the mechanism of geomagnetic activity by switching 

— * 

on the electric field E , after a period of time the velocity, v , and Joule heating, Q, 
will reach the values given by Eq. (21) and Eq. (23). 

(B) Hartmann Solution 

Cole [4 ] studied the problem of heating and dynamics near auroral electro- 
jets by considering the following model as shown below: 



4 

> z 


—4 


— * 


E ^ 


B 




— 1 i r~> 

-L L X 


The directions of the electric and magnetic fields are shown in the figure and the direction 
of flow is perpendicular to x - z plane. If the viscous effect of the fluid is taken 
into consideration, the momentum equation can be written as: 







( 24 ) 
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where j ~ o. (£ + v B ) . 
J x 1 y z' 

Equation (24) can be rewritten as 


A 


3 x 


b 2 v 


a, - 2 

— W v 


y 


Op 1 
T1 


E| | B 


(25) 


Let us consider the thickness of the sheet of the auroral zone to be 2L and 
assume that the motion of the gas outside this sheet is smallic We can then set the 
boundary condition as v * 0 at x = ± L. Then the solution of Eq. (25) that sati- 

y 

fies this boundary condition is given by 


where M ^ 
as follows: 


Ip I / cosh (M. x/L)n 
v (x) = fii ( 1 - 1 


B 


cosh M 


(26) 


j^| L (cf^ ) ^ • The heating within the sheet can be expressed 



cosh (Mj x/L ) 
cosh Mj 


(27) 


II -3 . Time Dependent Problem 

In the last section we have discussed steady solutions in magneto~gas-dynamics 
which indeed can help to understand the basic mechanism of Joule heating during a 
geomagnetic storm. However, they are not sufficient to explain the time variation of 
the geomagnetic disturbance. Actually all the quantities, electric field, magnetic 
field and velocity /are time dependent during the entire period of a geomagnetic storm. 
Therefore, the analysis for the unsteady case is necessary. 


Let us recall the equation of motion, namely Eq. (16) 


D v _ 
Dt 


=-V p +x+jxB 


Suppose that the current data j is known (by either observation or theory), then 
the induced magnetic field may be calculated by Eq. (8), together with the earth 


- 12 - 



magnetic field. Therefore, the last term of Eq. (16) becomes known. The problem 
is then reduced to a dynamic problem of the atmosphere with an additional given 
driving force (Lorentz force) and heat source (Joule heating). 

On the other hand, if the magnetic field data in the thermosphere are available, 
we will simplify the problem by eliminating the variables E and j from the system 
of equations. This is done by solving Eq. (18), namely, 

— ♦ 

£ = 4- - v x B +8 (]xB) . (28) 

2 -» a ] 

where a = a a ^ | B | . 

Substitution of j from Eq. (12) then gives 


E -■ -2-5-1 - ; x 1+ a 


9 x 8 


x B 


(29) 


when combined with Eq. (11) 'gives 


5 B 

at 


+ V x 


v X B 


- 9x ( v x B) + Vx 


a 

7 ~ 


(Vx B)x B = 0 . (30) 


We notice that Eq. (30) includes only two variables B and v . If B (x, y, z, t) 
is assumed to be given,then the velocity vector can be solved from Eq. (30) with 
properly imposed boundary conditions. 

Equations (11-18) may be reduced in the following form 


a o 


+ V . (p V ) = 0 


p^7= -*p + 


(9 x B)x B 


U 


D $ 


(7x B) 2 

— T~ 

Hb CT 


(31) 

(32) 


P 


Dt 


+ v* (xvT) -v* (p v) + cp . 


(33) 



This set of equations together with Eq. (30) formed a complete description of the 
time dependent problem. As observed from these equations, the only electromagneti 
quantity explicitly involved in this system of equations (30-33) is the magnetic field 

—4 

B . If B is given, this set of equations is reduced to ordinary dynamic equations. 



CHAPTER III 


NUMERICAL EXAMPLES 
III— 1 • Statement of the Problem 

A simplified model considers the thermosphere (80 - 200 KM) as a neutral 
inviscid atmosphere with finite electric conductivity. Hall effect can be neglected 
as we have shown in Chapter II. Thus, the set of fundamental equations is as 
fol lows: 


and 


nr-v*. - 


i-~ + 7 • (pv) 


= 0 


Dv _ 


Dt 


1 — » 1 — * — » 

— V P + — — - ( j x B ) 

p K P J 


D T — * — * 

p C div (X VT)+ j • E* 


E 1 ~ E + v x B 


7 x £ = - 


SB 

St 


j * o ( E + v x B ) 


p - P R T 


(34) 

(35) 

(36) 

(37) 

(38) 

(39) 

(40) 


We further assume that the atmospheric motions are confined to the vertical 
(z - axis) and horizontal (x-axis, E-W) directions. For the magnetic field, since 
observational data are notavailabl eat this stage, we assume quite arbitrarily that 

8 (x, y, t) = b (t) exp (- 8 r ) , 
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where r = \/ x + y and g is an arbitrary constant, which is the measure of 
the damping rate of geomagnetic variations. This assumed expression for the dis- 
turbed geomagnetic variations is such that the maximum disturbance is at the center 
of the region where we are interested, and it will decay exponentially according to 
the distance from the center, which also implies that the geomagnetic disturbances 
are confined in a finite region. Furthermore, the altitude dependence is also 
ignored in the present study. The quantity b (t) is taken from ground station magnetic 
data. In this calculation the May 1967 geomagnetic storm data recorded at College, 
Alaska; Dallas, Texas; and Honolulu, Hawaii were used. The 1964 Jacchia model 
atmosphere was used to calculate the pre-storm conditions of the atmosphere. The 
set of equations (34-41) was integrated numerically by using a finite difference tech- 
nique. The numerical process is explained in Diagram (2). It includes the following 
steps. 

— » 

1. Calculate E from Eq. (38) and (41). 

-* I 

2. Obtain E from Eq. (37). 

— » 

3. Calculate j from Eq. (39). 

4. Integrate Eq. (35), determine v for next time step, v (At). 

5. Integrate Eq. (36), calculate T (At). 

6. Calculate p (At) from Eq. (40). 

7. Determine p (At) from Eq. (34). 

8. Use v (At) to calculate E (At) from Eq. (37). 

Ill— 2. Results 

Numerical results are presented for winds, joule heating and electric field at 
the center of the storm and an altitude of 140 KM for May 1967. The calculations 
are based on the observational data recorded at College, Alaska; Dallas, Texas; 
and Honolulu, Hawaii. The Alaskan results show that the geomagnetic storm can 
generate the horizontal wind (East-West direction) (Fig. 4) on the order of maxi- 
mum 1000 M. sec. and vertical wind (Fig. 5) approximately an order of magni- 
tude smaller, which agree with the recent electrojet observation Fees [ 13 ]. The 

-3 -1 

joule heating (Fig. 6) is about a few erg -M - sec. and the electric field (Fig. 3) 
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on the order of 100 mV - M \ which agree well with the results given by Cole [ 14] 

and Wu, Matsushita/ and De Vries [15]. However, the calculated temperature in the 

present model is rather high. We believe this is due to the fact that the present 

model has not taken into account the effects due to viscosity. Considerably smaller 

winds, joule heating, and electric field result when the Texas and the Hawaiian data 

are used instead of the College, Alaska data. The horizontal (East-West) windv over 

-1 x 

Dallas, Texas (Fig. 10) is calculated to be . — . 1 00 M. sec. , and the vertical wind v 

_1 _2 z 

(Fig II) is~!0M. sec • The joule heating (Fig. 12) is approximately 10 erg - 

-3-1 - 1 

M - sec. and the electric field (Fig. 9) is on the order of 10 mV - M . The 

Hawaiian results show that v ~ lOM-sec. ^ (Fig. 16), v ~M-sec. ^ (Fig. 17), 

-4 -3 X _1 z _l 

Q ~ 10 erg - M - sec. (Fig. 18) and E~1 M- sec. (Fig. 15). The results 

show that the strongest thermal and dynamic responses in the thermosphere due to the 

magnetic storm are in the polar region and they become weaker as latitude decreases. 
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CHAPTER IV 


CONCLUDING REMARKS AND RECOMMENDATIONS 

In this study, a theoretical model for the dynamical responses to the geomagnetic 

storm in the thermosphere is derived from the magnetohydrodynamic theory. The 

validity of this proposed model is based on the assumption that the thermosphere behaves 

as an electric conducting fluid. This is true, because the degree of ionization in the 

-4 

thermosphere is ~ 10 and thus, the dynamical properties in the thermosphere are domi- 
nated by the neutral gas. In order to test the proposed model. A numerical example 
is presented. In this calculation, we have neglected the viscous effects and fixed our 
attention only on a narrow region of the thermosphere, and constant transport properties 
are assumed in the model. In general, the results we obtained are in good agreement 
with the wind measurements from an Agena satellite [ 13] and joule heating obtained 
by Cole [3 ] . However, the calculated density and temperature in this model are 
unreasonably high. We believe this is because we have ignored the viscous effects, 
the altitude dependence, and the gravitational wave effects. Therefore, we shall 
recommend that this model be improved in the following manner: 

1) By including the viscous effect. 

2) By including a global calculation model using a suitable spherical coordinate 
system. Thus, the latitude dependence can be incorporated into the model. 

3) Currently, our results are obtained for a point. We hope to do a calcula- 
tion to include whole region of the atmosphere from 90 Km - 250 Km, the 
structure of the winds and electric field can be obtained. 

4) By examining joule heating, the effects on gravitational wave. 

These will form a basis of our future studies. 
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